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We introduce matrix generalizations of the Navier-Stokes (NS) equation for fluid flow, and the 
Kardar-Parisi-Zhang (KPZ) equation for interface growth. The underlying fleld, velocity for the 
NS equation, or the height in the case of KPZ, is promoted to a matrix that transforms as the 
adjoint representation of SU{N). Perturbative expansions simplify in the TV — > oo limit, dominated 
by planar graphs. We provide the results of a one-loop analysis, but have not succeeded in flnding 
the full solution of the theory in this limit. 



I. INTRODUCTION 



■ The occurrence of complex probability distributions under non-equilibrium conditions is a challenging topic. Re- 
fSJ ' cently, a great deal of effort has been focussed on the study of the Kardar-Parisi-Zhang (KPZ) equation 

^. ^^^^V^h+^{Vhf+^M, (1) 

a deceptively simple prototype of non-equilibrium dynamics. The evolution of the field /i(x, t), denoting, for example, 
the height of a growing interface has a stochastic component governed by the noise ?7(x, t), usually assumed to be 
C ' Gaussian distributed, with zero mean, and correlations 

o 

Vn ; (r;(x, t)r](x', i')) = 2D8%^ - ii')S{t - t'). (2) 

ON 

The dynamical equation can be regarded as representing an elaborate filter that converts the simple correlations of 
the noise rj to complex correlations of the field h{x,t). For example, the two point correlations satisfy the dynamic 
scaling form. 



([/.(x, t) ^ /.(x', t')f) = |x - xf X/ (|x - x'l V|i - t'\) , (3) 



^ \ where z and x ^re the so called dynamic and roughness exponents. 

For A = 0, the linear diffusion equation gives z = 2 and x = (2 — d)/2. Equation (|l|) is rendered interesting and non- 
trivial due to the nonlinearity. The renormalization group (RG) flow of the effective coupling constant = X^D/v^, 
is given in a one loop perturbative calculation [hM] as 
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where iCd = Sd/i^irY, and Sd is the d-dimensional solid angle. Equation (||) suggests the presence of strong couphng 
behavior in all dimensions d. In d = 1 the dynamics is super-diffusive, characterized by exponents z = 3/2 and 
X = 1/2 JH]. At d = 2, the nonlinearity is marginally relevant, signaling strong coupling behavior; recent numerical 
studies [^indicate z w 1.612, and x ~ 0.386. In dimensions d > 2, a critical strength of A separates diffusive [z — 2 
and X = (2 — d)/2], and non-trivial regimes. Unfortunately, the strong coupling fixed point is not accessible by 
perturbative RG HJ^], which has recently been extended to two loop order (iQl- In an alternative approach, the one 
loop perturbation equations are converted into a set of self-consistent, so-called mode-coupling equations |p|-pl]|. The 
numerical solution of these equations remains controversial, with glassy solutions recently suggested as a possibility 
p^ . The mode-coupling equations have been shown to be exact when the field h is generalized to an N component 
vector, in a large N limit. 

The KPZ equation bears certain similarities (but also profound differences) with the Navier-Stokes (NS) equation 
for fluid flow. The fluid velocity field v(x, <) evolves according to 

^ + A(vV)v = -^ + i.V2v + f(x,t), (5) 
ot p 
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where v is the viscosity, p is the density, and A (=1) is introduced for book-keeping purposes. The pressure p(x, t) is 
adjusted to enforce the incompressibiUty condition V • v = 0. The qualitative difference with the KPZ equation is in 
the stirring force f (x, t), which has non-zero Fourier components f(k, uS) only at small fc ^ l/L, where L is the length 
scale at which energy is pumped into the system. The nonlinearity then transfers the energy into modes with higher 
wave-numbers A;, eventually dissipating it at the dissipative length scale. This establishes the Kolmogorov energy 
cascade in the intermediate (inertial) regime of wave- vectors A basic question in turbulence is the nature of this 
energy cascade, usually described by a scaling form for the energy spectrum, E{k) cx fc^^. (If the velocity correlations 
satisfy a dynamic scaling form as in Eq.(|^), then 

E = dkE{k) - J d'^kduo J d'^xdte'^^-'^-'^^'^ (v(x, i)v(0, 0)) 



and thus 



E{k) ~ k"^-^ j d'^xe'^ ^'x^^ 



and thus, as is well known, we have the exponenet relation C = 2% + 1.) The famous argument by Kolmogorov 
gives ^ = 5/3, in reasonable agreement with experiment. 

To apply the methods of dynamical RG it is usually assumed that f(k, cj) is Gaussian distributed with zero 

mean, and correlations 

(/i(k,L^)/,(k',L^')) = 2D{k)P,,{k){27rf+^d''{k + k')6{u; + iu'). (6) 

The transverse projection operator Pij(k) = 6ij — kikj/k"^ is a consequence of the condition V • f = 0. Thermal 
fluctuations in the fluid can be modelled by D{k) oc /c^, and lead to the familiar long-time tails in the velocity 
correlation functions More drastic forms of stirring are modelled by D{k) oc k~^ with —p < 2 In the problem 
of turbulence, the stirring at the longest scales resembles D{k) oc S'^{'k), which has the same scaling as p = d. This 
observation motivates some of the more recent applications of RG to turbulence jist . An extensive review of the 
application of RG to turbulence in given by Mou and Weichman in Ref. In particular, these authors develop a 
set of self-consistent equations for the problem, exact in an appropriate N oo limit, for which C = 3/2. 

The study of large N matrix theories started with the classic work of Wigner jl^ and has been developed [|l8|jl9[ 
over the years. Notable advances include the formulation of large N quantum chromodynamics pc|], and applications 
to random surfaces j2^. Intensive studies of the subject have continued in recent years |^^. Here we formulate 
and study large N matrix generalizations of the KPZ and NS equations by taking advantage of a crucial difference 
between vector and matrix models: the product of two vectors is not a vector while the product of two matrices 
is a matrix. Thus, we are able to promote the fields /i(x, t) and v(x, t) to N x N hermitian matrices and preserve 
the non-linear structure of the A'^ = 1 equations. As large N matrix models are generally quite different from their 
vector counterparts, we may hope to obtain results that are distinct from the earlier mode-coupling equations [ p^ . 
In particular, it is possible that techniques developed in large N matrix theory, such as the fact that its perturbative 
expansion is dominated by planar diagrams, may be brought to bear on this problem, leading to new analytic results. 
Although the primary interest is in the N = 1 limit, analytical insights may shed light on such controversial issues 
as the existence or absence of an upper critical dimension for the KPZ equation. The main focus of this paper is a 
perturbative analysis of the generalized KPZ equation introduced in Sec.|^. A brief analysis of the matrix NS equation 
is presented next in Sec. III. Prospects for further developments are discussed in the final Section, 



II. MATRIX KPZ EQUATION 

We promote the height /i(x, t) in Eq.(|l|) to a hermitian matrix field /i^(x, t), transforming as the adjoint represen- 
tation under the SU{N) symmetry group. The indices a and /3, which we call color indices, run over 1, • • • , A'^, with 
N large. Time evolution of the matrix is governed by the generalized KPZ equation, 

dh" A A A 

= jy.V^hl + v^S^pV^Kl^ + Y^h^Vh} + MVh'^pVKl^ ^5^V/i^, V/i^' + ^<S|V/i^V/i^; -f r;^(x, t). (7) 

All repeated indices are summed over; the contractions guarantee the SU (N) invariance. The stochastic noise is 
assumed to have zero mean, with correlations 

(r7^(x, t),7f;(x', t')) = 2 [d,52,5^; + D,5"X] '^'(^ " ^O'^Ct ~ t'). (8) 
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Although constrained by SU{N) invariance, the equation still admits 8 parameters' four at the linear level (i^i,2, ^1,2) 
and four nonlinearities (Ai^2,3.4)- Due to this large number of parameters, Eqs.(|3) are quite formidable. It would 
be convenient if we could focus on a subset of parameters such as (z^i, Ai,Di). To check if such a simplification is 
possible, we appeal to a graphical perturbation expansion. The diagrammatic approach |^ to the KPZ equation may 
be generalized to the matrix represented by two points: The propagator is a double line that keeps track of the 
flow of the color indices a, /?, as well as momenta and frequencies. The nonlinearities are represented as triplets of 
double lines, while averaging over noise joins pairs of such lines. These diagrammatic entities are depicted in Fig. (|^). 
By substituting these entities for the corresponding elements, each term in the original perturbative expansion now 
has several counterparts in the matrix theory. Summing over the internal labels generates a factor of N for every 
closed loop. Keeping only the diagrams with the largest power of N will hopefully lead to some simplification. 

To gain insight into the structure of the perturbation series, we started by setting all parameters except {i>i, Ai, Di) 
to zero, and performing a one loop RG analysis. From the perturbative expansion of the propagator, and the correlation 
function, we obtain respectively 

diyi_, Kd{d-2) NXlDi dv2 _ , ^, ^Kd{d-2) XfU^ 

— -{Z- 2.)Ul — -2 — ^[Z- Z)U2 



dl ^ ' Ad vi dl ' U 
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dl A ul M 4 vl 

The equations for vi and Di are similar to those of the scalar KPZ, with 2A^ replaced by iVA^. We also see that V2 
and D2 are already generated at this order, although these corrections are smaller by a factor of 1/N . The reductions 
by powers of 1/N arise because to create the pairing of indices (the 'U' turns) corresponding to 1^2, etc., out of 
(1^1, Ai, Ui), the lines have to cross at some point. However, as is well known in matrix theory, such non-planar 
diagrams carry smaller powers of N, and can be neglected in the N ^ 00 limit [po| . The next question is whether 
these generated parameters, albeit small, can feed back into the recursion relations for (i/i,Ai,-Di) significantly, so 
that their inclusion is necessary. Clearly, replacing any of the original graphical elements with the above generated 
ones results in a smaller power of N from their smaller magnitude. The question is whether the rearrangements of 
the indices can generate compensating powers of N. For all simple graphs that we examined the answer was negative. 

We thus assume that the subspace (Ai = X/^/N, vi = v,Di = D) is closed up to order and study the matrix 
equation 

dh°' A 

^ = i.V2/i^ + ^V/i;;V/.^ + ry^(x,t), (10) 

with a noise term described by 



(,7^(x,i)r/f;(x',t')) - 2i?C<'5'(x-x')5(t-t')- (H) 



(We shall shortly demonstrate that this choice of scaling the parameters with N , which is in line with the one loop 
result, leads to a consistent large N theory at all orders.) The reduced matrix equation is invariant under the 
transformation 

A . . \ 2j 



h'^{^,t)^h^p (^x+-^ui,<j +5^u.x + 5^^u^i, (12) 

which generalizes the so called Galilean invariance |^ of the scalar KPZ equation. Referring to Eq.(^, we see that 
the dynamic and roughness exponents are defined by requiring that the effective equation preserves its form under 
the transformation x — > 6x, t —* b^t, and h b^h. Under these rescalings, and with 6 = e^, the non-linear coupling 
A evolves according to 

^ = (x + z-2)A, (13) 

quite generally for any N. To preserve the invariance in Eq.(^2|) under the above rescalings, the coefficient A must 
remain constant, leading to the exponent relation 

X + z = 2. (14) 

This identity holds at any fixed point with finite A. 
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In d = 1, Eq.(|T^) also satisfies a fluctuation-dissipation condition which implies that the probability distribution 



Volh] (X exp 
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(15) 



is a stationary solution of the Fokker-Planck equation for V{[h],t) In deriving the matrix version of the Fokker- 
Planck equation we must be careful about the non-commutativity of matrices. Going through the usual formal steps, 
we obtain 



d^x { tr— 
oh 



2VN 



Dtr 



Sh- 



(16) 



The nesting of the parentheses in the first term of this equation requires some explanation: the large round brackets 
contain all the matrix quantities to be traced over {V is of course not a matrix), while the square parenthesis indicates 
that S/6h should also act on V. Inserting Eq.(15), wc find that the right hand side of Eq.(|^) becomes, for a general 
dimension d, proportional to J (i'*xtr[(V/i)^V^/i]. For d = 1 the integrand is a total divergence, indicating that Vq 
is a stationary solution. For this stationary state x — 1/2 (and z — 3/2 from Eq.(|l^), also in agreement with the 
RG results of Eq.(|^). Thus, in d = 1 the exponents for = 1 and N ^ oo are identical; also a feature of the vector 
equations [Q. In this context, it is interesting to note a recent study of a trimer deposition model on a line that is 
claimed to be equivalent to an = 2 matrix model p5[ |. Numerical results on this model appear to suggest z sa 2.5. 

By integrating over the noise in a functional integral description, we can also obtain a more traditional field 
theoretical formulation of the problem as 



{Z) = 



vh / vtj j[h]s i - j/v^/i;^ 
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(17) 



where J[h] is the Jacobian associated with the transformation from h to 77. The value of J[h\ depends on the 
discretization of the evolution equation. In an Ito choice of discretization in which the noise 77 (x, t) only affects the 
field at ft,(x, t + t), the Jacobian associated with |9r/(x, t)/3/i(x, t + T)| is simply a constant |^^. We shall thus ignore 
this factor henceforth. The (5-functions in Eq.(pT|) can be implemented by using a conjugate field h as in Rcf. [Q. 
Integrating over this field, as well as the noise 77, then leads to the action 



S{h) 



1 
2D 



d'^xditr 



dh 

at 



A 



2ViV 



(18) 



We can now remove two of the remaining parameters by the rescalings, t t/v and h ND /v h, leading to 



Sih)^^ J d'^xdttr(^(^|-V2)/7-|(V/7)^ 



(19) 



with g as defined after Eq.(^). With its cubic and quartic interactions, the above action is reminiscent of that of a 
non-abelian gauge theory such as quantum electrodynamics without quarks. 

The Feynman rules corresponding to Eq.(p^) are given in Fig.(|^). Each propagator line has a factor of 1/A^, while 
each vertex is proportional to N, and each loop generates a factor of N. Thus, a given diagram with E external lines, 
/ internal lines, V3 three-point vertices, V4 four-point vertices, and L loops, is associated with a factor of with 
P = L + Vs + V4^ — I. Using the standard topological identity L = I — V3 — V4 + 1, we obtain P = 1. For E — 3 and 
E ^ A; this certainly coincides with the bare scalings of three and four point vertices in Eq. (p^) . Thus the choice of 
scaling the nonlinearity with ^/N indeed leads to a consistent large N expansion. 

We also carried out a one-loop perturbative RG, directly on the action in Eq.(p^. The basic idea is to integrate 
out short wavelength modes, with A/6 < fc < A where A is a cutoff. The effective action for the remaining modes 
then has a leading gradient expansion of the form 



Sih) 



N 

y 



(i''x dt tr 



h- 



(20) 



depending on two parameters a and /3. The same coefficient a multiplies both dth and the nonlinearity, as required 
by the Galilean symmetry in Eq.(^2|). To calculate the effective action, we only have to look at terms quadratic in h. 
Thus, it suffices to evaluate the renormalized propagator. 



4 



Note that a term proportional to tr(V/i)^ is also allowed by symmetries, and is in fact generated. However, such 
a term can be eliminated by transforming to a moving coordinate frame h h + ct. Thus, in calculating the one- 
particle- irreducible self-energy function E(fc,aj), we only have to extract the coefficients of the and k"^ terms in a 
low frequency and wave- number expansion. The one-loop corrections to the propagator are indicated in Fig. (^. In 
fact, the two diagrams are identical, except that one is smaller by a factor of due to the crossing of lines. We 
shall keep track of both diagrams by including a factor of (1 + l/N), thereby also correctly reproducing the = 1 
case A simplifying feature is that the quartic interaction term ~ g^iyh^ in Eq. (|l9|) does not enter into the 

calculation to this order. Evaluating the frequency and momentum dependence of the remaining self-energy diagram 
is cumbersome, but straightforward, and leads to 



I](fc,w) 
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Adding the self-energy to the bare propagator, we obtain the renormalization parameters, 
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lA/b {2ttY 



We have evaluated the integral in Eq. (|2|) over an infinitesimal shell by setting b — 1 + di and writing 

Kddi, with Kd = SdA'^~'^ /{2tt)'^, where Sd is the d-dimensional solid angle. 

The cutoff in the effective action can be restored to its original value by setting x = 6x', accompanied by t — h^t' 
and h = b^h' . The parameters in the quadratic part of the renormalized action can be reset to unity by choosing the 
exponents 
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Finally, the flow of the coupling constant is given by dg / d£ — (x + 2; — 2) g, resulting in 

,2 / 1 \ 11J2 
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„ , g- / 1 \ lld^ -5d- 18 
'^'+-6''^['+n) did + 2) 



(24) 
(25) 

(26) 



Surprisingly, this flow equation looks quite different from that obtained by standard dynamic RG methods 
which generalize Eq.(0) to 



dg^ 2 
-dl=' 



, 9 l\2d-3 



(27) 



The resolution to this discrepancy may be due to the fact that the flow equations are not fundamental. The two 
equations in fact agree for both d = 2 and d = 1. In d = 1, the RG is constrained to give the correct exponents due to 
the symmetries embodied by Eqs.(|l2|) and (|l5|). It is indeed easy to check that by substituting Ki{l + 1/N)g*^ = 4 
in Eqs.(p^ we recover z ~ 3/2 and x — 1/2. The exponents at the phase transition between the weak and strong 
coupling phases in d = 2 + e dimensions are also correctly given by z = 2 + O(e^), x = + O(e^). (The correlation 
length at this transition diverges with an exponent v = e^^ + 0(1).) Thus all ph ysi cal quantities calculated from 
Eqs.(p6|) and (27) appear to be identical. The coefficient of the quadratic term in Eq.(p6|) changes sign at d = 1.5265 as 
opposed to d = 3/2 for Eq. (|2^) . However, there is probably no physical significance to this dimension, as it disappears 
in the two loop calculation of Ref . 0| . Even simpler fiow equations (without any higher order terms) are obtained in 
an RG of the KPZ equation that proceeds from a mapping to directed polymers pi]. 



III. MATRIX NS EQUATION 

For the Navier-Stokes equation we promote the velocity vector field v — {vi, ■ ■ ■ , Vd) to a set of x hermitian 
matrices v^(x, t), imposing the incompressibility diVi = as before. The non-linear term has four matrix counterparts 
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Vj{djVi), {djVi)vj, Vj{diVj), and {diVj)vj. We will suppress the matrix indices henceforth. (We shall neglect various 
linear and nonlinear terms that can be constructed using combinations of the unit matrix and/or trv as in Eq.([7h. 
They are again expected to be irrelevant in the large N limit.) We further require the matrix generalization of Eq.(||) 
to satisfy two important physical conditions: (1) a suitable generalization of Galilean invariance, and (2) conservation 
of energy in the absence of dissipation. Subject to the above limitations, we arrive at 



dt 



{vj,djVi} + 77-7= [vj,djVi + diVj] 



2VN 



2VN 



diP 
P 



vd]v. 



(28) 



where the curly ({...}) and square ([...]) brackets indicate, respectively, matrix anti-commutation and commutation. 
(The matrix products ensure SU{N) invariance.) The commutator term proportional to k has no scalar counterpart. 
The 'pressure' matrix p{x.,t), is again adjusted to enforce the incompressibility condition diVi = 0. Finally, it is 
assumed that the random (matrix) force is Gaussian distributed with zero mean, and correlations that in Fourier 
space are given by 



(29) 



The transverse projection operator again follows from incompressibility, and difi = 0. 

A form of Galilean invariance is satisfied by the generalized equation, if in a frame of reference moving with velocity 
u, the generalized velocity transforms to 



(x,t) = v;^ I x + ^ut,t 



(30) 



It is easy to check that Eq.(E8p is invariant under this transformation. A natural choice for the generalized energy of 
the fluid is -E = p/2 J d'^xtr(ViVi). In the absence of dissipation (ly ~ 0) and forcing (f — 0), the change in energy is 
governed by 



dE 
'dt 



^ P J d'^xtr (vidtVi) = - j d'^xtr (^ ^^J— Vi {vj,djVi} + ^^ "^i K": c'jO'i + diV^\ + Vid^p^ . (31) 



Using the incompressibility condition djVj = 0, and reordering the matrices inside the trace, we can transform the 
above expression to 



— = - I d'^xdAx 
dt ' ^ 



Xp 



2VN 



ViViVj + Vjp 



-0, 



(32) 



i.e. energy is conserved as required. In what follows we will set k to zero for simplicity. It would be interesting to 
study how k flows under the renormalization group. 

The above conditions are sufficient to give the scaling exponents for Eq.(p8|) under most forcing conditions. Just 
as before, we define the exponents x ^-^d z through the correlation function 



(t,.(x, t)v,i^', t')) = <5., |x - xf ^/ (|x - xT /|< - t'\) 



(33) 



As in Eq.(|l4|), the requirement of generalized Galilean invariance leads to the exponent relation x~'r z — 1- We need 
one more exponent relation to determine x s^nd z separately. 

Let us first examine thermal noise with D{k) = Dk^ (model A in the language of Ref. j|]). In this case, the 
linearized equation with A = has a steady state probability distribution 



PoM oc exp 



2D 



(34) 



It is then straighforward to show that the contribution of the non-linear term (A 7^ 0) to the probability current of the 
Fokker-Planck equation is the integral of a divergence. The manipulations establishing this result are almost identical 
to those in Eqs.(|3l|) and (|32| ). Thus by ensuring energy conservation, we have also set up a fluctuation-dissipation 
condition j23|. From Eq.(|34|) we can read off the dimension of the velocity as x = — d/2, leading to z = 1 -I- d/2 from 
Galilean invariance. (Note that in d = 1, Eq.(p8|) for model A is simply the derivative of Eq.([ic|).) 

Now consider a more general stirring force with D(k) — Dpk~P. Because the nonlinearity is proportional to the 
gradient, the perturbative series generates only powers of fc^. Thus the coefficient Dp, the most relevant component 
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of noise as fc — > 0, is not renormalized to any order in perturbation theory. This non-renormahzation leads to 
the exponent relation X — {z — d + p)/2. Together with the constraint from Galilean invariance, we obtain 
X = {1 — d + p)/3 and z — (2 + d — p)/3, and a Kolmogorov exponent of C = 2% + 1 = (5 — c? + p)/3. For the 
case of large scale stirring {p — d), the Kolmogorov result of ^ — 5/3 is recovered. However, by appealing to this 
scaling analysis, we have simply reproduced previously known RG results, and gained nothing new from the matrix 
equation. By contrast, the large N equations introduced by Mou and Weichman fl^ can be summed exactly in a 
large N limit. These authors find that the perturbative exponents break down for —p < 1 — d, when the exponent z 
sticks to 1. The non-linear term now dominates dt'v, and equating the bare dimensions of v • Vv and the noise leads 
tox=(l^c? + p)/4. According to Ref. for p > d, this leads to the exact result of ^ = 3/2 for the Kolmogorov 
exponent in their large N generalization of the NS equation. Unfortunately, we have not yet succeeded in finding 
the exact behavior of our matrix generalization as iV ^ oo. Mou and Weichman do consider an adjoint SU{N) 
generalization, presumably equivalent to Eq.([2^). However, we believe that the construction of the equation, and its 
analysis, is more transparent in our matrix formulation. 



IV. DISCUSSION 



In this paper, we introduced and analyzed matrix generalizations of the KPZ and NS equations. It is clearly possible 
to construct similar generalizations of other dynamical equations, such as for the time dependent Landau-Ginzburg 
process. Our original motivation was to find a set of closed form expressions that are exact in the large N limit, 
and which hopefully yield different results from those based on other N oo generalizations of these equations 
p0| , pl]p6| . The large N limit does indeed lead to certain simplifications of perturbative expansions; most notably in 
the dominance of planar diagrams However, we have yet to succeed in taking advantage of this fact to reduce 
the perturbation series to a closed set of self-consistent equations. Nevertheless, we are hopeful that progress towards 
this goal can be made. A number of techniques for dealing with large TV matrix theories have been developed over 
the years. For instance, consider the model where a single N x N hermitian matrix ip is taken randomly from the 
probability ensemble P{ip) oc exp[—NtTV{ip)]. While difficult to evaluate diagrammatically, quite remarkably, this 
model can be solved for any V, by using the orthogonal polynomial approach It is an intriguing possibility 

that some analogous approach may be developed for the type of problems discussed here, and for large N quantum 
chromodynamics (QCD). Another approach towards solving large N QCD involves the concept of a master field. In 
the present context, this amounts to finding a master field /imaster which dominates the functional integral in Eq.(pT|). 
Certainly, the KPZ problem is considerably simpler than QCD, since we have one master field instead of the four 
associated with A^, = t,x,y,z. However, we do not have gauge invariance, which is powerful enough to constrain 
the master fields in QCD to be independent of space-time coordinates. Yet another approach is to consider an RG 
flow ||3^ in TV, attempting to relate the theory for TV -I- 1 to the one for N. The hope is that the model will flow 
towards a fixed point as ^ oo. In any case, the field of large N matrix theories is evolving rapidly with recent 



advances in pure mathematics [^l 32 1. Some of these advances may ultimately prove fruitful in our context. 

The calculations in this paper were limited to one loop perturbation theory. Surprisingly, in the case of the KPZ 
equation, distinct flow equations for the effective coupling constant were obtained from two different starting points 
(the dynamic equations, and an effective field theory). However, the two RG schemes are equivalent at the critical 
dimension, and appear to result in the same exponents. We did not carry out a similar program for the NS equation, 
as simple scaling arguments appeared to be sufficient for obtaining exponents. It would, nonetheless, be interes ting 
to check if a direct analysis reveals any surprises in this case. By the same procedure as in Sec.||, we can turn Eq.(|2q) 
into a field theory. The complicating new feature is the presence of the new interaction term proportional to k. We 
don't know of any symmetry condition that will keep k = under renormalization. We also note that the interaction 
terms, while still quartic in the velocity fields, now contain only two derivatives, rather than the four in the KPZ case. 
It will interesting to study the RG flows in the (A, k) parameter space. The previous scaling analysis will be invalid 
at fixed points where either parameter is infinite (or both are zero). 

Finally, we conclude by noting yet another mapping of the KPZ equation; to directed polymers in random media 



|p3| . For a scalar field /i(x, t), the transformation C/(x, t) — exp [ghix, t)/2], changes the action in Eq.(|l9|) to 

S{U)^— j d'^^dtir{U-^dtU ~U-^dluf . (35) 

Promoting U to an A x A^ matrix leads to another possible generalization. We note that this is not, strictly speaking, 
equivalent to Eq.(^|) due to the non-commutativity between matrices. However, we can take this expression as 
another generalization, no less valid than Eq.(p^, of the KPZ equation. The above action resembles a non-relativistic 
version of non-linear sigma models, and principal chiral field models, discussed in the field theory literature, with the 
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crucial difference that the matrix U is not unitary. Perhaps some of the methods developed in this area can also be 
brought to bear on the problem at hand. 
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FIG. 1. The elements of a diagrammatic expansion starting with the dynamical equation. 
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FIG. 2. The elements of a diagrammatic expansion starting with the field theory action. 



FIG. 3. The diagrams contributing to the one loop correction of the propagator in the action. 



9 



